An interesting charmonium state formation and decay: pp — > 1 D 2 — > 1 P± r y 
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ABSTRACT 

Massless perturbative QCD forbids, at leading order, the exclusive annihi- 
lation of proton-antiproton into some charmonium states, which, however, have 
been observed in the pp channel, indicating the significance of higher order and 
non perturbative effects in the few GeV energy region. The most well known cases 
are those of the 1 5'o (rj c ) and the 1 P±. The case of the 1 _D 2 is considered here 
and a way of detecting such a state through its typical angular distribution in 
the radiative decay 1 D 2 — > 1 Pi'j is suggested. Estimates of the branching ratio 
BR( 1 D 2 -> pp), as given by a quark-diquark model of the nucleon, mass correc- 
tions and an instanton induced process are presented. 
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Several hadronic two-body decays of charmonium states are forbidden within 
leading order perturbative QCD [1,2], but are nevertheless observed to occur with 
decay rates comparable to or even bigger than those of allowed decays [3]; the 
most well known examples are the J/ip — > VP [4], and the 1] C — > BB, VV [5] 
channels, where P is a psdeudoscalar meson, V a vector meson and B a baryon. 
Indeed the observed decay rates for J/ip — > p-zr, K*K and r\ c — > pp, pp, 00, K*K* 
are difficult to explain within conventional perturbative QCD. Recently, also the 
1 Pi coupling to pp has been established [6] , despite being equally forbidden by the 
helicity conservation rule of massless pQCD [7] . 

Among the attempts to solve these problems, non leading contributions [8], 
two quark correlations inside baryons [9], quark mass effects [10] and gluonic con- 
tents of the J/ip [4] and the p c [11] have been considered. Higher order Fock states 
might help with the J/tjj — > pir decay [8], but their contributions to other pro- 
cesses are not clear; diquarks and mass corrections do not help much with the rj c 
forbidden decays, whereas gluonic contributions seem to be more promising [12]. 
Recently a dynamical model for such contributions, with instanton-induced non 
perturbative chiral symmetry breaking, has been used to obtain a good agreement 
with the data on r(p c — > pp) [13]. 

We consider here yet another case of a charmonium decay which should be 
forbidden according to perturbative QCD, namely 1 D 2 — > pp. Its observation 
would be very interesting because, among the above non perturbative mechanisms 
invoked to explain the other forbidden decays, no one seems to be able to account 
for a sizeable decay rate: as we shall see, both mass corrections and diquarks give 
very tiny decay rate values and instanton induced processes are strongly suppressed 
with increasing Q 2 values [13]. 

Let us briefly recall why the 1 D 2 — > pp decay is forbidden in massless pQCD. 
This charmonium state has quantum numbers J =2 ^: parity, angular mo- 
mentum and charge conjugation conservation only allow a final pp state with or- 
bital angular momentum L = 2 and spin S = 0. S = implies that the p and the 
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p must have, in the charmonium rest frame, the same helicity, which is forbidden 
by the QCD vector coupling of hard gluons to massless quarks and antiquarks. 
Such helicity selection rule can only be broken by terms proportional to m q /m c 
or hx/nic, where m q and m c are respectively the light quark and charmed quark 
masses and kr is the quark intrinsic transverse momentum. The current masses 
of u and d quarks are very small compared to the charmed quark mass and terms 
proportional to m q /m c are indeed negligible; terms proportional to /ct/w c might 
be more relevant, but no comprehensive treatment of these contributions, together 
with other higher twist effects, has yet been performed. 

Let us now consider the 1 D 2 state created in pp annihilations, choosing the 
z-axis as the proton direction in the pp centre of mass frame. It is then clear, from 
what we said above, that the 1 D 2 state can only be created with the spin third 
component J z = 0; such charmonium state is then purely polarized and its spin 
density matrix has only one non zero component, 



This peculiar property reflects into the decay angular distributions of the 1 D 2 . 
One radiative decay which is expected to be observed with a large branching ratio 
is 



which is dominated by an electric dipole transition. The angular distribution of 
the photon, as it emerges in the rest frame of the 1 D 2 , is then simply given by [14] 



where 9 is the photon polar angle and an integration has been performed over the 
azimuthal angle. 

The observation of such an angular distribution in pp exclusive annihilations 
should then be a clear signal of the formation and decay of the 1 D 2 state; the full 



p 00 ( 1 D 2 ) = l. 



(1) 




(2) 



W 7 (0) = -(5 -3 cos 2 9), 
8 



(3) 
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chain of processes to be looked for, according to the observed or expected decays 
of the 1 Pi state [6], is: 

pp -> X D 2 -> x Pi7 -> (??c7) 7 -»■ (777) 7 , (4) 

or 

-> 1 D 2 -> X Pi7 -> ( J/^tt°) 7 -> (e+e-yr ) 7 . (5) 

Notice that the expected mass of the 1 D 2 state is M D = (3788 ± 7) MeV [15]. 

So far the 1 So(r) c ), 3 Si(J/i^ and J 7 ), 3 Pi (Xci) an d 3 P 2 (Xc2) charmonium 
states have been observed to couple to pp; the corresponding branching ratios, 
BR( 2S+1 Lj -> pp), are typically of the order of 10 -4 to 10 -3 [3]. Curiously, 
the rj c — > pp branching ratio, which should be zero according to lowest order 
perturbative QCD, is among the largest ones. Recently, also the 1 Pi has been 
observed in the pp 1 Pi — > J/tfj7r° channel [6], with an estimate for the product of 
the two branching ratios BRQ-P 1 -> pp) BR( 1 P 1 -> J/^tt°) ~ 10" 7 . Notice that, 
similarly to what explained for the 1 _D 2 , also the x Pi decay into pp is forbidden 
by leading order pQCD [7]. The 3 Po state has not yet been observed, but this is 
presumably due to its small (< 10 -2 ) branching ratio into J/^7; this makes the 
full process through which one looks for such a state, pp — > Xco — ► J/^l e + e~7, 
a difficult one to detect. The analogous situation for the 3 Pi and S P 2 states is 
much more favourable in that their branching ratios into J/i/j'f are respectively 
~ 0.27 and 0.13 [3]. 

Thus, it is natural to expect a 1 D 2 — > pp branching ratio similar to that 
observed for the other charmonium states. However, this would be very difficult 
to explain; to see why, we now briefly consider several possible non perturbative 
contributions to such process. 

Mass corrections to 'forbidden' charmonium decays have been considered in 
Ref. [10] for rj c , Xco — > PP', they yield sizeable values of T(x c o — > pp), but very 
small ones for r(?y c — > pp), actually a factor ~ 10~ 4 smaller than data. Following 



the same procedure and notations as in Ref. [10] we have computed the helicity 
amplitudes for the decay x Di — > pp; the only non zero ones are: 



A ++ , M (e) = -a. 



M9) = § ] fl^R>\0)§ F e(l-4e*f 



d 2 

a M 



pi pi— X2 pi pi — 2/2 

((9) / dx 2 \ dx 3 / dy 2 / dy 3 
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(6) 



[x 2 y2 + (x2 - 2/2) 2 e 2 ] [1 + x 2 y2 -x 2 -y2 + {x 2 - ?/2) 2 e 2 ] 

1 1 

[x 3 y 3 + (x 3 - y 3 ) 2 e 2 } [(1 - x 2 )y 3 + (1 - x 2 - y 3 ) 2 e 2 ] 

: 73 (X2 - ?/2) 3 X 

[X2V2 ~ 2 {x2 + 2/2) + (x 2 - y 2 ) 2 e 2 \ 

V? x (231)v? y (321) - ^(132) [^(321) + p„(123)] 

- [^(132) + ^(231)] ^(123)] (I-X2- y 3 ) 
^ x (123)^(213) - ^(321) [^(213) + ^(312)] 

- [^(321) + <^(123)]<^(312)] (1 - x 2 ) 
+ [^(213)^(123) - ip x (312) fo,(321) + p„(123)] 

- [^(213) + <^ (312)] <^ (321)] (1 - x 2 ) 

where M D is the 1 D 2 mass and R"(0) is the value at the origin of its wave func- 
tion second derivative. <p z (i,j,k) = <p(zi, Zj, Zk) denotes the proton distribution 
amplitude and we refer to Ref. [10] for further details. Here we only notice that 
e is the ratio of the proton to the charmonium mass, e = m p /M D , so that in the 
massless limit, e — > 0, indeed A±±-m = 0, as required by perturbative QCD. 
From Eq. (6) one obtains the decay rate 

:i-4e 2 ) 1 /2 



T( 1 D 2 -> pp) = 



E 



d(cosO)\A Xp , XpM m 



= ^ 7 r^|^(0)| 2 |F Ar |V(l-4e 2 ) 9 / 2 ^ 



(7) 
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where I is the multiple integral appearing in Eq. (6). 

From Ref. [9]* also the decay rate of the 1 D 2 into two gluons can be obtained: 

rCD 2 ^gg) = ^ §r\R"(0)\\ (8) 

D 

By assuming the total hadronic decay rate of the 1 D 2 to be approximately given 
by Eq. (8), one can get an estimate of the branching ratio BR( 1 D 2 — > pp) by 
taking the ratio of Eqs. (7) and (8), so that the unknown quantity R"(0) cancels 
out. The result strongly depends on the choice of the distribution amplitudes 
<f(xi, X2, xs); according to the different choices adopted in Ref. [10] one obtains 

BRCD 2 -> pp) ~ 10" 8 -T- 10" 12 . (9) 

Eq. (9) clearly shows how mass corrections could not account for the eventual 
observation of the 1 D 2 — > pp decay; the small values obtained for the branching 
ratio are mainly due to the factor (x 2 — y 2 ) 3 contained in the decay amplitude, 
Eq. (6). This is similar to what happens for the 1] C — > pp process, where mass 
corrections are also very small, due to a factor (x 2 — y 2 ) in the amplitude [10]; in 
the present case, actually, the situation is even worse, because of the third power 
of (x 2 — y 2 ). In fact, in the r\ c case, mass corrections lead to BR{r\ c — > pp) ~ 
10~ 6 4- 10~ 10 [10], a result far away from the observed value BR{r\ c — > pp) ~ 10~ 3 , 
but bigger than the values given in Eq. (9). 

One can similarly show that also two quark correlations could not explain a 
branching ratio for the 1 D 2 —> pp decay of the order of 10 -4 ; a vector diquark 
component of the proton allows the decay, by allowing helicity flips at the gluon- 
vector diquark coupling [9], but, once more, the numerical values turn out to be too 
small. This can be explicitly checked by repeating the same procedure followed 
above for mass corrections; the expression of the decay helicity amplitudes, in 
the quark-diquark model of the proton, can be found in Ref. [9] and, again, it 

* Notice that in Ref. [9] the 1 D 2 state is named f 2 
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contains a small factor (x — y) 3 . One finds, with little dependence on the choice 
of the distribution amplitudes, 

brCd 2 -> pp) ~ 10" 8 . (10) 

Among other non perturbative effects proposed to explain unexpectedly large 
branching ratios, the presence of the fundamental (L = 0) trigluonium states, 
with quantum numbers J = h , 1 , 3 , in the 3 GeV mass region, has 
been proposed [4,11]. The first two states, mixing respectively with the r\ c and 
the J/ip, might explain some of their 'mysterious' decays. However, a similar 
explanation for the 1 _D 2 , the presence of a 2~+ glueball with a mass close to 3.8 
GeV, looks much less natural and realistic. 

Let us consider finally the instanton induced mechanism proposed in Ref. [13] 
for the r\ c — > pp decay: we know that its contribution decreases very rapidly with 
increasing Q 2 and, indeed, already for the decay of rj' c , with a mass ~ 3.6 GeV, is 
much smaller than for the rj c [13]. Considering the still higher mass of the 1 _D 2 , 
M D ~ 3.8 GeV, we cannot expect this non perturbative contribution to be large 
enough to produce a branching ratio for the process X D 2 — > pp similar to those 
observed for the other charmonium states. 

We have thus seen how several possible non perturbative effects cannot con- 
tribute significantly to the 1 D 2 coupling to pp; on the other hand we know that 
leading order pQCD predicts no coupling at all, whereas higher order corrections 
are difficult to evaluate and have never been computed. A similar situation oc- 
curs with the 77c, with the difference that for such particle one might expect a 
significant gluonic contribution [12,13]. Therefore, the eventual observation of a 
BR( 1 D 2 — > pp) ~ 10~ 4 , analogous to the values observed for all other charmonium 
states which can couple to pp, would pose an intriguing challenge to the theory. 

The 1 D 2 state could be looked for in the mass region M D ~ 3788 MeV [15] 
and in the reactions suggested by Eqs. (4) and (5), which should exhibit a typical 
decay angular distribution of the 7 in the first step of the process. In fact the 1 D 2 
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created in pp annihilation is in a pure spin state with J z = and its decay into 
1 Pi7, dominated by an El transition, has the simple angular distribution given 
in Eq. (3). Hopefully, such radiative decay has a large branching ratio, so that 
the processes of Eqs. (4) and/or (5) can be detected. This is not unrealistic if 
one notices that the 1 L>2 state, due to its expected mass and quantum numbers, 
cannot decay into pairs of D and/or D* mesons. 
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